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Abstract 

The coherent scattering of photon in the Coulomb field (the Delbriick 
scattering) is considered for the momentum transfer A <C m in the frame 
of the quasiclassical operator method. In high-energy region this process 
occurs over rather long distance. The process amplitude is calculated tak- 
ing into account the multiple scattering of particles of the intermediate 
electron-positron pair in a medium. The result is the suppression of the 
process. Limiting cases of weak and strong effects of the multiple scattering 
are analyzed. The approach used is the generalization of the method de- 
veloped by authors for consideration of the Landau-Pomeranchuk-Migdal 
effect. 
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1 Introduction 



The nonlinear effects of QED are due to the interaction of a photon with electron- 
positron field. These processes are the photon-photon scattering, the coherent 
photon scattering (the elastic scattering of a photon by the static Coulomb field 
called often the Delbriick scattering), the photon splitting into two photons, and 
the coalescence of two photons into photon in the Coulomb field. Among these 
processes the coherent scattering of photon has been observed and investigated 
experimentally (see reviews flU and @) during last decades, the most accurate 
measurement of the coherent scattering has been performed not long ago in the 
Budker Institute of Nuclear Physics 0. The photon splitting into two photons 
in the Coulomb field was observed for the first time at Budker INP only recently 
||]. Observation of photon-photon scattering is still a challenge. 

History of the coherent photon scattering study can be found in the men- 
tioned reviews. There is a special interest to the process for heavy elements 
because contributions of higher orders of Za (Z\e\ is the charge of nucleus, 
e 2 = a = 1/137, h — c — 1) into the amplitude of photon scattering are very im- 
portant. This means that one needs the theory which is exact with respect to the 
parameter Za. The amplitudes of the coherent photon scattering valid for any 
Za for high energy photon (u ^> m) and small scattering angle (or small momen- 
tum transfer A) were calculated in H, Q. The approximate method of summing 
of the set of Feynman diagrams with an arbitrary number of photons exchanged 
with the Coulomb source was used. Another representation of these amplitudes 
(in the Coulomb field) was found in 0, using the quasiclassical Green function 
of the Dirac equation in the Coulomb field. This Green function in a spherically 
symmetrical external field was obtained in M where the coherent photon scat- 
tering in the screened Coulomb potential was investigated as well. Lately this 
Green function was calculated for "localized potential", and the coherent photon 
scattering was analyzed using it 0. Recently the process of the coherent photon 
scattering was considered [IT| in frame of the quasiclassical operator method (see 
e.g. |L(J) which appears to be very adequate for consideration of this problem. 

The coherent photon scattering belong to the class of electromagnetic pro- 
cesses which in high-energy region occurs over rather long distance, known as the 
formation length. Among other processes there are the bremsstrahlung and the 
pair creation by a photon. If anything happens to an electron while traveling this 
distance, the emission can be disrupted. Landau and Pomeranchuk Jl2| showed 
that if the formation length of bremsstrahlung becomes comparable to the dis- 
tance over which a mean angle of the multiple scattering becomes comparable 
with a characteristic angle of radiation, the bremsstrahlung will be suppressed. 
Migdal |HJ developed a quantitative theory of this phenomenon which is known 
as the Landau- Pomeranchuk-Migdal (LPM) effect. 



A very successful series of experiments (see |M, |L5]]) was performed at SLAC 



during last years. In these experiments the cross section of bremsstrahlung of 
soft photons with energy from 200 KeV to 500 MeV from electrons with energy 
8 GeV and 25 GeV is measured with an accuracy of the order of a few percent. 
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The LPM was observed and investigated. These experiments were the challenge 
for the theory since in all the previous papers calculations are performed to loga- 
rithmic accuracy which is not enough for description of the new experiment. The 
contribution of the Coulomb corrections (at least for heavy elements) is larger 
then experimental errors and these corrections should be taken into account. 
Recently authors developed the new approach to the theory of the LPM 



effect in frame of the quasiclassical operator method. In it the cross section 
of bremsstrahlung process in the photon energies region where the influence of 
the LPM effect is very strong was calculated with term oc 1/L , where L is 
characteristic logarithm of the problem, and with the Coulomb corrections taken 
into account. In the photon energy region, where the LPM effect is "turned off', 
the obtained cross section gives the exact Bethe-Heitler cross section (within 
power accuracy) with the Coulomb corrections. This important feature was 
absent in the previous calculations. The contribution of an inelastic scattering of 
a projectile on atomic electrons is also included. We have analyzed (see |p~8|] ) 



the soft part of spectrum, including all the accompanying effects: the boundary 
photon emission, the multiphoton radiation and the influence of the polarization 
of a medium. Perfect agreement of the theory and SLAC data was achieved 
in the whole interval of measured photon energies. Very recently we apply this 
approach to the process of pair creation by a photon in [fJU] . 

In the quasiclassical approximation the amplitude M of the coherent photon 
scattering is described by diagram where the electron-positron pair is created 
by the initial photon with 4-momentum k\ (cu,ki) and then annihilate into the 
final photon with 4-momentum h.2- For high energy photon u 3> m this process 
occurs over a rather long distance, known as the time of life of the virtual state 

where q c > m is the characteristic transverse momentum of the process, the 
system h — c — 1 is used. When the virtual electron (or positron) is moving in 
a medium it scatters on atoms. The mean square of momentum transfer to the 
electron from a medium on the distance If is 



2 4 

q 2 = A7iZ 2 a 2 n a Llf, L(q c ) = ln^, q 2 mm = oT 2 + A 2 + ^, (1.2) 



where a = e 2 = 1/137, Z is the charge of nucleus, n a is the number density of 
atoms in the medium, A is the photon momentum transfer (A = |k 2 — ki| <C q c ), 
a is the screening radius of atom. 

The coherent photon scattering amplitude M can be obtained from general 
formulas for probabilities of electromagnetic processes in the frame of the quasi- 
classical operator method (see e.g. [TO). It can be estimated as 



au q 2 a 2 



M ~^T— ^ = — <£■ (1-3) 
2irl f n a q£ irn a 
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We use the normalization condition ImM = uoa p for the case A = 0, where a p is 
the total cross section of pair creation by a photon. 

In the case of small momentum transfer q s = ^Jq 2 Cmwe have in the region 
of small A <C m 

q 2 = m 2 , M —In -. 1.4 
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At an ultrahigh energy it is possible that q s ^> m. In this case the char- 
acteristic momentum transfer q c is defined by the momentum transfer q s . The 
self-consistency condition is 

2 2 2huj Z 2 a 2 n a L(q c ) 
q c =q s = - 2 , (1.5) 

where L(q c ) is defined in (|1.2|). So using ( |1.3[ ) one gets for the estimate of the 
coherent photon scattering amplitude M (the influence of the multiple scattering 
manifests itself at the high photon energies such that m 2 a/uj <C 1) 



2Z 2 a 3 u I A? 



M-^-yin-^^, A^ = y/2nuZ*a*n a » A 2 . (1.6) 

In the present paper we use consider the influence of the multiple scattering 
on the process of the coherent photon scattering for A -C q c . The theory of the 
coherent photon scattering in the Coulomb field in frame of the quasiclassical 
operator method is stated in Sec. 2. We give here an alternative presentation of 



approach developed in [IT] in a more formal way. In Sec. 3 we apply the method 
developed in [^6[ , |T7[| to investigation of influence of the multiple scattering on 
the process of the coherent photon scattering. The general formulas describing 
this influence were derived and the asymptotic cases of the strong and weak 
effects are analyzed. 



2 Coherent scattering of a photon 
in the Coulomb field 

2.1 Formulation of approach 

The coherent scattering of a photon in the Coulomb field (the photon with mo- 
mentum ki = (a>i, kx) — > k 2 = {uJ2i ^2)) is represented by the electron loop in the 
Coulomb field, and uj\ = u>2 = oj. The corresponding amplitude is 

T — — % J d 4 xi J rf 4 x 2 Tr [ei exp(iA; 1 xi)G'(xi, x 2 )e2 exp(iA; 2 x 2 )G'(x2, Xi)] , 

(2.1) 
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where e = e^7 M , e M is the photon polarization vector, G{x\,x 2 ) is the electron 
Green function in the Coulomb field the standard representation of which is 




G(x 2 , Xl ) = { " n w } ^(-) ' • (2-2) 

Here ^^(xi) are the solution of Dirac equation in the Coulomb field, signs (+) 
and (— ) relate respectively to positive and negative frequencies. 

As well known, see e.g. [[21], Sec. 12, in the noncovariant perturbation theory, 



which we use here, in high-energy region (lo ^> m) the diagram contributes into 
the amplitude T where the electron-positron pair is first created by the initial 
photon with 4-momentum k\ and polarization e\ in time t\ and then annihilate 
in time t 2 > £i into final photon with 4-momentum k 2 and polarization e 2 , while 
the contribution of the interval in which t 2 < t\ is of the order m 2 /uj 2 . Such 
contribution is neglected in the scope of our method. Taking this into account 
and substituting the Green function into the amplitude ( |2.1|) we find 

T = i £ / dti / dt 2 §{t 2 - ti)V nm (ei, k ls ti)V* m (e 2 , k 2 , t 2 ), (2.3) 

where 

V nm (e,k,t) = J d 3 r^ + (x)aeex P (-ikx)¥~\x), (2.4) 
It is evident that V nm (e, k, t) ( j2.4|) is the matrix element of pair creation by a 



photon in an external field. In the quasiclassical operator method (see [[10], Sec. 3) 
the wave functions in an external field can be presented in the form 

¥±\r,t) =< r\¥ ± \~P,H)exp[Ti(H±eip)t} \n >, H = Vm 2 + P 2 , (2.5) 

where <p = <p(r) is the potential of an atom, \n > is the state in the configuration 
space at the time t — 0, the functions ^^(P, Ti) have form of wave functions for 
free particles in the momentum space (# (±) (p,e). We substitute (^TsJ) into 



and take into account completeness of states |r > (/ |r >< r\d\ = I). After this 



we convey consistently in ( p.4[ ) the operator exp(zkr) to the right up to \m >, 
and then the operator exp [i(Ti + e<p)t) from < n\ to exp [i(Ti(k — P) — e(p)t). As 
a result we obtain the combination of operators 



V nm (e,-k,t) = (n\^ + (P(t))(xe¥-\k-P(t)) 

x exp(i(Ti - ecp)t) exp(i(H(k - P) + &p)t) exp(ikr) \m) (2.6) 



the disentanglement of which can be performed in the standard way (see (T(J . 
Eqs.(3.7)-(3.11)): 

L p (t) = exp [i(Ti - &p)t\ exp [i(H(k - P) + e<p)t] 



exp(iut)T exp 



. } kP(t) ' 



o 



, P=(Ti,P), (2.7) 
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where T is the operator of the chronological product, and the expression 



R p = # (+) (P, n)e^-\k - P, H(k - P)) 



(2.8) 



becomes the Heisenberg operator depending on time P = P(t), Ji — H(t). It 
can be expressed in terms of two-component spinors. Substituting these results 
into Eq.fl2.4p we find 



V„ m (e,k,t) = (n\R p (t)Texp 



kP{t) 



LU - H(t) 



dt 



Rp(t) = (A(t) - iaB(t)) <f s , A 



2H(lo - H) 



exp(ikr)|my , 
e(kxP)), 



B 



2H(uj - H) 



[emuj + (eP)(k-2P)] 



(2.9) 



where tp$, (p s are two-component spinors describing polarization of created positron 
and electron. Here the Coulomb gauge is used. 

Substituting ( |2.9|) into (|2.3| ) and performing summation over states \n > and 
\m >, as well as summation over spin states s, s we have the amplitude of the 
coherent scattering of a photon in the Coulomb field 



T 



2na 



S 2 i = Tr 



ijdhj dt 2 $(t 2 -t 1 )S 21 
m t{^) exp[— iAr]m 1 (t 1 ) 



A = k 2 - ki, 



(2.10) 



where the operations Tr means sum of the diagonal matrix elements both in the 
configuration and the spin spaces, 



m lj2 (t) = m(e 1>2 , k x>2 , t) = (A(t) - ier~B(t)) Texp 



ki, 2 P(t) 
u - Kit) 



dt 



(2.11) 



functions A(t), B(t) are defined in ( [2.9|) . Let us remind that P(t) is the momen- 
tum operator in the Heisenberg picture: 



P(i) = exp{-iHt)Pexp(iHt), H = H + V, 



(2.12) 



where V = eip 

We present the evolution operator as 



exp(-ifft) = exp(-iHt)N(t), N(t) = exp(iHt) exp(-i(H + V)t) (2.13) 
Differentiating the last expression over time we find 



= -i exp(i?#)V(r) exp(-z(W + V)t) = -iV(r + vt)N(t), v = ^ 
dt ri 



(2.14) 
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The solution of this differential equation for the initial condition N(0) = 1 is 



N(t) = Texp 



-i J V(r + vt')dt' 



(2.15) 



If the formation length of photon scattering is much longer than the charac- 
teristic length of electron scattering (u/m 2 3> a), one can present the dependence 
on time of the operator P(t) in fl2.9|) as (see e.g. fllQ| , Sec. 7.1) 



P(t) = (-t)P(-oo) + tf(t)P(oo), P(±oo) = iV+(±oo)P(t)iV(±oo). (2.16) 
Let us introduce states convenient for further calculations 



\i >= iV + (-oo)|p, >, |/ >= iV+(oo)| P/ >, 
P|Pi >= Pi I Pi >, Pp/ >- P/ P/ > 



(2.17) 



The states \i > and \ f > are the eigenvectors of the operators P(— oo) and P(oo) 
correspondingly 

P(-oo)|i >= iV + (-oo)PN(-oo)iV + (-oo)|pi >= pi\i > 

P(oo)|/ >= iV+(oo)PiV(oo)iV + (oo)|p / >= p f \f > . (2.18) 



The general expression for the amplitude of photon scattering (Eqs. (|2.10 ) 
and (|2.11|) ) is simplified significantly for small momentum transfers A = |A| 
(A<m, A/e < I/7). With the accuracy ~ A/m we can put that the operator 
P commutes with exp[— iAr] and substitute the vector k2 in the expression for 
1^2(^2) by the vector ki. Using ( |2.16| ) and the states \i > and \f > ( |2.18| ) for 
construction of the matrix element we have 



/|m(P(t))exp 



--Ar 

2 



™ fi (e,p(t)) (f\exp 



-Ar 

2 



m fi (e, p(t)) = [^(-t)m(e, p { ) + i?(t)m(e, P/ )] , 



(2.19) 



where we describe motion of particles of the created pair as a trajectory in "the 
form of an angle" . Bearing in mind that the commutator 



(2.20) 



we can discard operator T in the expression for operator N(t) ( p,15|) and with 
relativistic accuracy (i.e. with accuracy up to terms ~ m/uS) present the operator 
N(t) as 



N(t) = exp 



J V(g,z + t')dt' 



(2.21) 



where the axis z is directed along the momentum of the initial photon ki, q is 
two-dimensional vector transverse to axis z. 



7 



Using ( 2.17[ ) an d Q2.21 ) we calculate the matrix element in ( [2. 19 ) 



/|exp 



-lAr 

2 



J < f\r >< r\i > exp 



--Ar 

2 



d 3 r 



YJ ~ ¥) 6 ( q|1 " 2) ' 



(2.22) 



where 



< /|r >=< p / |AT(oo)|r >= iV(oo) < p/|r >, < p/|r >= 



exp(— ipfr) 



< r\i >=< r|iV + (-oo)|pi >= N + (-oo) < r|pj >, < r|p; > 



(2vr) 3 / 2 ' 
_ exp(zp^r) 

(27r) 3 /2 ' 



1 

/(Q) = — / exp[iQg + i X {Q)]d 2 Q, x(Q) = ~ I V(q, z)dz. (2.23) 
lixi J J 



We introduced here vectors 



Pi = p, p/ = p + q, 



(2.24) 



the function /(Q) is the scattering amplitude in the eikonal approximation. 
Substituting Eq. fl2.19| ) into T 2i Eq.( |2.10| ) we have 



S 2 , = Tr 



m 



(t 2 ) exp[-iAr]mi(ti) =Xl Tr m jU e2 ' P(*2))w/t(e 1 , p(ti)) 



i,f 



x ( ?| exp 

_ tf(A||) 
" (2vr) 2 



*Ar 

2 



I/) (/I exp 



--Ar 

2 



d 3 p Tr mt(e 2 ,p(t 2 ))^/i(ei,p(ti)) 



x / / f <u - ^r) /* f q± + J , ?n = o. 



(2.25) 



2.2 Helicity amplitudes 

In the above analysis we traced the transition to the expressions calculated on the 
trajectory of particle "in the form of an angle" in the momentum space. Precisely 
these trajectories determine the amplitude of the coherent photon scattering for 
uj/m 2 > a (see Eq.( pTT0|) ): 



1 a 



2tt5(A||)) J^j-J da{q u A)R 21 , (2.26) 



-21 



(2tt)* 



where 



dff(qj_, A) = / L - ^) f* U± + ^ ) rfqj 
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R 2 i = dh dt 2 '&(t2 - ti)Tr m^(e 2 , p(t 2 ))m/i(ei, p(*i)) 



1 

2/^1/ ^2i?(*2-<i)>C(e2,ei;v2±,vij.)exp -i^ ^ kiv{t)dt 



UJ 



£(e2,ei; v-ll, vu.) = 



(e^m^Xeinvix) + — (e^ei 



(a; - 2s) s 



(e2V2±)(e!V 1± ^ 



Vi, 2 ± = V±(t 



1,2 



p ^1 

u = - = (l,v), e' = u-e, n = — , (vi j2 ±n) = 0, 

e a; 



(2.27) 



where (env) = (e(n x v)). 

We calculate now R21 in this expression for the trajectory in "the form of an 
angle" (see (gig)): 

p(t)=pi?(-t) + (p + q)#)- (2-28) 

We get for this case 
m 2 ^ = (e*ng)( ei ng) + ^^|^-(e*g)(e lg ) + (e* ei )^; 



Pj 



m + (p± + qj 



2 • 



m + pi 



#0 



m 2 + (pj 



q± 



m 

-t 5 ^.2.29 

) 2 m 2 + p^ v 7 



The functions g and go have very important properties: 

g(q± = 0) = 0, O (q± = 0) = 0. 



(2.30) 



In the case of complete screening when the screening radius a <C u/m 2 the 
amplitude of photon scattering T21 ([2.26D is imaginary at any A. Taking into 
account that at A = this amplitude is connected due to unitary relation 



with the known probability of pair photoproduction, see e.g. pT|, p2| , we will 
calculate the difference 

«?2i = T 21 (A)-T 21 (0). 

The interval | q_i_ | < | A | <C m contributes into this difference and one can expand 
the functions g, go as a power series in q±: 



g 



q± 



2p±(q±P_ 



m + Pj_ [m z + pjj 



2m(q ± p_L) 
(m 2 + p^) 2 ' 



(2.31) 



We introduce the dimensionless variable u = p±/m and carry on averaging over 
angles of u (the azimuthal angle of the component of electron momentum in the 
plane which is perpendicular to the direction of initial photon n) using formulas 



M 2 

v~v~ = — UiUjU k ui = — (dijdki + 5 ik 5ji + 5u5 jk ) 
2 o 



(2.32) 
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where Sij is the two-dimensional Kronecker delta. Substituting Eq.( 2.3l|) into 
( P^D and using ( ggg ) we get 

m C 4 #2i = (e2nq ± )(e x nq ± ) + — (^qj.)(eiqj.) 



+q 2 L (e 2 ; e 1 ) 



1 



2ee r 



(C 5 



4 c- 1 

or 



, C = 1 + u 2 . (2.33) 



It is convenient to describe the process of photon scattering in terms of helicity 
amplitudes. We choose the polarization vectors with helicity A 



e A = -j= (ex + i\e 2 ) , e 1 = v 



e 2 = n x v, A = ±1, 



e x e* x = 1, e x e*_ x = 0, e A x n = i\e x , (2.34) 
where n is defined in Eq.( [2.27 ). There are two independent helicity amplitudes: 

M ++ = M__, M + _ = M_ + , 

where the first subscript is the helicity of the initial photon and the second is 
the helicity of the final photon. When the initial photons are unpolarized the 
differential cross section of scattering summed over final photons polarization 
contains the combination 



2[|M ++ | 2 + |M 4 
Substituting (|2.34| ) into ( |2.33| ) we find 

2ee' 



(2.35) 



m 4 C 4 



i Ik 2 - i) • — — rK i) 



2ql ee> 

. — -008 2^ cos^=-. 



We define the amplitude of the coherent photon scattering M 2 i as 



7T 



T 21 = -5(A||)M 21 , 



a; 



then the cross section of the coherent photon scattering has a form 



da. 



27t5(A|| 



\T, 



21 



, d 3 k 2 



167T 2 



|M 2 i| 2 (ifi 2 



^2.36) 



^2.37) 



^2.38) 



Substituting ( p37| ), (gggfl and ( gg§ ) into Eq.(ggg) we obtain for the helicity 
amplitudes 



u 00 



M xx , = l ^-m 2 J dejd(j da(q ± ,A)R xy 



[2.39) 



where AA' is ++ or H — . For the chosen normalization of the helicity amplitudes 
the unitary relation with total cross section of pair photoproduction a p is 



-Im M 4 



(7(7 • 
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e e 



CT p . 



^2.40) 



So we present the helicity amplitude M ++ as 

M ++ = 5M ++ + iua p , 

w oo 

5M ++ = % -^-m 2 J de J d( J (da(q_ ± , A) - dcr(q ± , 0))R++. (2.41) 
o 1 

Since the forward scattering amplitude (A = 0) with helicity flip vanishes we 
have for amplitude M + _ 



U) oo 



6M + _ = ^-m 2 JdeJdcJ da(q ± , A)f2+_. (2.42) 



o 1 



Integrals over qj_ in Eqs.( [2.4lD and fl2.42| ) we denote 5q 2 ± and q^cos2<y9. We 
calculate these integrals using the scattering amplitude /(Q) in the eikonal ap- 
proximation Eq.( [2.23D . We find 



5q 2 =Zx = J qi(d<7(qj_, A) - da(q ± , 0)) 
= / rfqxql f* U± + 5f\ f (<U - ^ ) - /* (<li) / (<Li) 

oo 

= 2tt| ( X '(f?)) 2 (Jo(Af ? ) - l)f?df?, 
o 

oo 

-q 2 cos2^ = Z 2 = 2n J (g)) 2 J 2 (Ag)gdg. (2.43) 

o 

For the screened Coulomb potential 

U( r ) = — exp(-r/a), x'(^) = ^^-Ki f- 
r a Va 



Z, -!.ZVi;(^), F a (z)= 2 ^±l f ln(z + VTT^) 1. 
2 / ^vl + z z 



Z 2 = 4irZ 2 a 2 F 1 ( — ) , F x (z) = 1 



Aa 



L= In (z + VTT^ 2 ) .(2.44) 

_|_ ~2 \ / 



The functions Fi(x) and i^(^) encounter in the radiation theory and in the theory 
of Landau-Pomeranchuk-Migdal (LPM) effect. In this derivation the following 
relations have been used 

oo oo i / t / n \ \ 

v 7 



1 2 



Kq{x)K x (x)xJ x ((3x) 
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1 1 00 J 1 00 

2 + £ / ^oOr)) 2 *Ji(/*r) = - - J (K (x)) 2 J ((3x), 



oo 

y xKl{x)xdx = -. 



(2.45) 



The cross section of e e + pair photoproduction in the case of complete screen- 
ing (a <C u/m 2 ) to within terms ~ m/u has a form (see e.g. Eq.(19.17) in pip 



9m 2 



ln(ma) + - - /(Za) - — 



where 



/(0 = Re[^(l + 20-^(l)]=e 2 E 



(2.46) 



(2.47) 



^ n(n 2 + £ 2 ) ' 

here ip(x) is the logarithmic derivative of the gamma function. Substituting 
(gg|)-(ggg) into Eqs. p^IQ and ( ggg) we obtain 

l l 
da; y fixx'fxx', 



where 



Im M xx , 



4Z 2 a 3 u 



(2.48) 



o o 



= 1 — 2x{l — x) + 4x(l — x)y{l — y), = x(l — x)y 2 

/++ = (^) + ln(ma) + ~ - /(Za) - i 

= ln (™) - ^=f ln (* + - + §, 

Here we passed to the variables x — e/ui, y = 



(2.49) 



2.3 The coherent photon scattering in different cases 



The important property of Eq.( |2.48| ) is that the dependence on the screening 
radius a originates in it from the Born approximation. In this approximation in 
the case of arbitrary screening the radius a enters only in the combination 



+ ?||> 911 



m 2 uj 



2ee' 



m 



x{\ — x)y 



a? (2 ' 50) 



Since in the course of derivation of Eq.( |2.48| ) we did not integrate by parts over 
the variables x and y, we can extend Eq. fl2.49Q on the case of arbitrary screening 
making the substitution 



1 

a 



q 2 + a~ 2 



Qef, 



2q ef 



(2.51) 
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Recently the amplitudes of the coherent photon scattering were derived in [ I I 
for arbitrary interrelation between parameters in the Moliere potential. 

We calculate the real part of scattering amplitude using the dispersion re- 
lations method. We can apply this method directly since in the region of mo- 
mentum transfer under consideration (q± <C m) the dependence of amplitude 
T 2 i oc M 21 /uj on q\\ has the form (see Eqs.(CT)-(ETl5l), Q), O) 

T{q\\) = Jdtxj dhFituhWh-tJexpi-mih-h)}. (2.52) 



Thus the real and imaginary part of the amplitude M 2 \/uj are connected by the 
Sokhotsky-Plemelj formulas 



F f 1m 
Re / (z) = — I —j 

71 J Z 1 



(2.53) 



Using this transformation and Eqs.( |2.48|) , ( |2.49| ) we obtain 



where 



1 1 



4Z 2 a 3 u , 

Re M\\i = — — - — j dx dy /ixx'Vxx', 







(2.54) 



7T 



/aa'(^ + ?||) ~ fx\'(\z - q\\\) 
fxx'(q\\(z + l))-fxx'(q\\\z-l\) 



dz 

z 

dz 

z 



(2.55) 



We consider now different limiting cases depending on interrelation between 
a and uu/m 2 . In the case of complete screening (a uo/m 2 ) we can use directly 
Eqs.( p.48|) , ( p.49| ). In this case the functions f\y are independent of x and y and 
the corresponding integrals are 



1 1 



1 1 



J dx J dy = 7 -, J dx J dy /z + _ = 





For the scattering amplitudes we have 

28Z 2 ct 3 a; 







Im M 4 

Re Mi 



9m 2 
0, Re M. 



Im M 4 




2Z 2 a 3 m f 
9m 2 



(2.56) 



(2.57) 



The results obtained are consistent with Eqs.(33) and (36) of || where calculation 
has been done for the Moliere potential. 
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In the case a ^> uj/m 2 (the screening radius is very large, or in other words 
we consider the photon scattering in the Coulomb field) we have, taking into 
account Eqs. (^49|) , fl2iHD 



i m 

/++ = In — • 

911 Wl + s 2 



2S ' + 1 lui.s, + w'l i - 



1 



In ( s c + x/1 + s 2 cl . 



41 

/(^) + -, 

= — = —x(l-x)y. 2.58 
2gy wr 



The expressions for amplitudes M ++ and M H in this case are consistent with 

the results obtained in ||. At A = we find known result (see Eq.(8.1) in ||) 



.A 



In 



■2lo 
. m 

oo 



-x(l — x)y 



1 41 
- - - HZ a) + — , 

2 v ; 42 



7T 



- In 



z + 1 dz 

\z — 1| -2 



- / In 

7T J 




1 + Z(i2 7T 

= 2' 



.28Z 2 a d u; 

i 

9m 2 



1 — z z 
109 



m M ; 42 



7T 



M+_ = 0. (2.59) 



At A = and arbitrary interrelation between a and uj/m 2 we get 



f 1 m 1 

A + = ln----/(^) + -, 



1 

27 



In 



<5 2 (z + l) 2 + q; 2 (l-x)y tfe 
<5 2 (z — l) 2 + x 2 (l — x) 2 y 2 z 



am 



m a 

20J \ r 



[2.m) 



The photon scattering amplitude in this case for arbitrary value of parameter 5 
was found recently in |§. 



3 Influence of the multiple scattering on the co- 
herent scattering of a photon 

3.1 Basic equations 

When a photon is propagating in a medium it dissociates with probability oc a 
into an electron-positron pair. The virtual electron and positron interact with a 
medium and can scatter on atoms. In this scattering the electron and positron 
interaction with the Coulomb field in the course of the coherent scattering of pho- 
ton is involved also. There is a direct analogue with the Landau-Pomeranchuk- 
Migdal (LPM) effect: the influence of the multiple scattering on process of the 
bremsstrahlung and pair creation by a photon in a medium at high energy |T^| . 
[ |j~3|j . However there is the difference: in the LPM effect the particles of electron- 
positron pair created by a photon are on the mass shell while in the process of 
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the coherent scattering of photon this particles are off the mass shell, but in the 
high energy region (this is the only region where the influence of the multiple 
scattering is pronounced) the shift from the mass shell is relatively small. To 
include this scattering into consideration the amplitude R21 Eq.( |2.26| ) should be 
averaged over all possible trajectories of electron and positron with the weight 
function <icr(q, A). This operation can be performed with the aid of the distri- 
bution function averaged over the atomic positions of scatterers in the medium. 
This procedure was worked out in details in fl7j (Sec. 2), (Sec. 2). The am- 
plitude of the coherent photon scattering for a photon propagating in a medium 
can be derived in the same way as Eqs.(2.4)-(2.6) of JT7J: 



M(e 1; e 2 ,A) = -g- J d 3 pexp (-i~r) J dV j dr'C (e 1; e 2 ;tf',tf) 



(2tt)%, 
xF (r', v', r; r, v) exp 



i- ki(r' — r) 



(3.1) 



The distribution function F (r', v', r; r, v) satisfies the kinetic equation 
dF fr', v', t) 



,dF(r'y,r) 



dr dr' 
n a J da(v', v", A) [F (r' ; v", r) - F (r', v', r)] 



(3.2) 



and the initial condition 

F (r, v, 0; r, v) = 5(r - r')5(v - v), 

here n a is the number density of atoms in a medium and in the case of the 
screened Coulomb potential 



da(v',v",A) = d<x(q,A), q = £ (v'-v) 
The combinations entering in Eq. (|2.27| ) are 



(3.3) 



/ kv(t)dt = lut- ki(r(t 2 ) - r(ti)) -> lut - k x (r' 

ti 

£ (ei, e 2 ; v 2 , vi) -»• £ (ei, e 2 ; t?', •&) 



(3.4) 



We can proceed with further calculation (integration over v', r') of ( |3 . 1| ) 
by analogy with the procedure used in Eqs.(2.7)-(2.18) of |17| |. As a result we 
obtain for the differential cross section of the coherent photon scattering with 
the multiple scattering taken into account 



^^iMKe.A) 



(3.5) 
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where 



M e 1; e 2 , A = / — / dtexp -it e* e] >o (0,t) 

n a J ee J 

o L 

- t ( ei nV)(e;ny(0, t)) - i ^ ~ ^ (ei V)(e^(0, t)) 



(3.6) 



In derivation we have changed variables into t, £? as in Eq.(2.7) of |Tj|. The 
function (p^(Q,t), tp^ = (<p ,ip) satisfies the equation 

■° r ~ 1 ' -Hipp, H = p 2 - i V(Q,A), p = -iV e , 



dt 

V(q, A) = f (1 - exp(zq^)) / L 



±) r L + £ U, (3.7) 



with the initial conditions 

yo(e, o) = 5(e), o) = p5(g). 



(3.8) 



The potential V(g, 0) = V(q) was used in the theory of the LPM effect (see 
(Sec.2)): 



N ^i( T ,4 nr \ ^ 2-nZ 2 a 2 ee'n a _ , a; 2 
= Qf? 2 Ui + ln— -2C) , Q= — ^, L 1 = ln^ 



m 4 c<j 



^ = lSSZ-^e-/, / = / (Z«) = (Za) 2 E 



1 



A 



J fc(fc 2 + (Za) 2 ) 



(3.9) 



where C = 0.577216 ... is Euler's constant. We can restrict ourselves to calcula- 
tion of the difference of the potentials only: 

AV(q,A) = V{q,A)-V(q) 

= {- \ (1 - exp(zqe)) (d(7(q, A) - da(q, 0))d 2 q, (3.10) 

turn 4 J 



where the cross section dcr(q, A) is defined in Eq.( [2.26| ). In the above formulas g 
is the two-dimensional space of the impact parameters measured in the Compton 
wavelength A c which is conjugate to the two-dimensional space of the transverse 
momentum transfers q measured in the electron mass m. 



3.2 Amplitudes of coherent photon scattering under in- 
fluence of the multiple scattering 

When parameter A<m the main contribution into the integral in ( |3.10| ) gives 
the region |q| < A/m <C 1. For g < 1 one can expand the integrand in ( |3.10|) in 
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powers of q£. Using Eq.( 2.23 ) we get 

ee'n a r (xlq) 2 



AV(Q, A) = ^ / ^(x'(x)) 2 (exp(-^x) - l)d\ 
cum 4 J x 1 



x( x ) 



A, 



Za 



exp — 



x 
r 



2ZaK n x 



X = X 



111Z" 



-1/3 _ ^s2 



exp / 



m 



(3.11) 



where K (z) is the modified Bessel function. Integrating AV over azimuthal 
angle of the vector x we obtain 

oo 

AV(g, A) = -2Qq 2 J Kf(x) (1 - J {px) + J 2 {/3x) cos2y?)xcfe 



2Qg 2 / K 2 (x) cos 2ip + Kf(x) (1 - J (f3x)) xdx 



Qq 2 



F 1 (^\ cos2^ + F 2 f^ 



, (3 = Aa s 



(3.12) 



where <p is the angle between the vectors q and A, J n (z) is the Bessel function, 
the functions F\{z) and F 2 {z) are defined in Eq.( |2.44 ). Some details of calculation 
are given after Eq. (|2.44|) . 



(3.13) 



Summing AV(g, A) (gig ) and V(t?) (gj) we get 
V (q, A) = Qg 2 f Li + In ^ - 2C - Fx f cos 2^ - F 2 f ^ 



Starting from Eq.( [3.6|) and advancing as in Sec. II (Eqs.(2.2)-(2.12)) of |17 
we find 



M. 



2iam 2 uj f de 



++ 



— f% (o| ai (g; 1 - Go- 1 ) + , 2P (g; 1 - Go 1 ) p|o 



M. 



^ / £ (0|* (elp) (G- - G-) (e + p) |o) , (3.14) 



+- 



where 

2 — I — £ /2 2 pp' 

si = 1, s 2 = 5—, s 3 = ^, G S = W + 1, G = p 2 + 1, (3.15) 

7i is defined in Eq. ( |3.7| ). In derivation it is convenient to use properties of eA 
given in ( |2.34|) . Note, that coefficients si,s 2 coincide with the same coefficients 
in the theory of the LPM effect for pair creation by a photon [J20| . As in 
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(Eq.(2.17)) and @ (Eq.(2.5)) we split the potential V(g, A) as 

V{q, A) = V c (g, A) + v(g, A), K(e,A) = g p f? 2 , q p = QL s , 

v( e , A) = (2C + In ^ + Fi ^ cos 2^ (3.16) 

In accordance with such division of the potential we will write the photon scat- 
tering amplitudes in the form 

M++ = M c ++ + AfQ, M+_ = M%_ + Mj 1 !, (3.17) 

where M+ + is the contribution of the potential V c (g, A) and is the first 

correction to the scattering amplitude due to the potential v(g, A) (see |16| 



Eqs.(2.28) and (2.33), and |g Eqs.(2.10) and (2.15)). Since dependence of the 
potential V(g, A) on £ is the same as in [^(J and the expression for M ++ (|3.14| ) 



formally coincides (up to external factor) with Eq.(2.2) in [20] the result for M ++ 
can be taken from |20| (Eqs.(2.10) and (2.15)) with substitution q — > q p , L c — > 
L s . The combination of photon polarization vectors entering in the amplitude 
M+_ ( |3.14| ) gives in evaluation the expression (e*_'V^j (e + g) = 0. Substituting 

this combination in calculation of the first correction M+l (see Eqs.(2.31) and 
(2.32) of |T(J) we obtain expression of type ( e _f?) i e +Q) — g 2 exp(2iip) . So in 
integration over the azimuthal angle (p only the term with cos2<£> survives in 
formula for v (g, A). As a result we get for the main term 

am 2 u r de , , s MC 



am uj f de , , . 

m c ++ = - — / —M v »)> 

27rn a J ee 



$ s (iO = v s I dte 



oo 

1 1\ / 1 1 



-it 



\ sinh z z) 2 V sinh 2 z z 2 



= Si (\np-ip (p+ ^)) + s 2 ^V(P) -lnp+ ^ , (3.18) 

where v s = 2yfiq~ s , z = u s t, p = ij (2z/ s ), i]){x) is the logarithmic derivative of the 
gamma function. 

The first corrections to the amplitudes are defined by 

++ Aun a L s J 66' sy sh + ~ Aim a L s 1 \2 J J 66' 3V sh 

dze . z/ s (rfs 9\ z )\ 



o 

/i(z) = ^ln £> 2 + In — — In sinh z — C^j g(z) — 2 cosh zG(z) 
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F 3 (u a) 



tu s s 3 



dze 



( \ 53 



1 + ^-pC(2,p) 



smb. z 2 

o 

g(z) = z cosh z — sinh z, t = t\ + t2, z — u s t, p 

z 

G ( z ) = /(! - ycothy)dy 
o 

= -f-^- ta ( 1 - e " 2 ") + ^( e ""). 



z 

2zZ 



(3.19) 



here Li 2 (x) is the Euler dilogarithm, ((s, a) is the generalized Riemann zeta 
function. Use of the given representations of functions F 3 (v s ) and G(z) simplifies 
the numerical calculation. 

In the region of the weak effect of scattering (| v s | <C 1 , g c — 1) the interval 
z <C 1 contributes into the integrals ( |3.18[ ), ( |3.19| ). In this region 



F *(y) 



Sl I 6 + ~60~ + 126 



+ S2 y + lT + -63" 



S 2 - Si 



F«(v) 



v 



1 + \v 2 + ^ 4 ) s 3 . (3.20) 



Substituting these expressions into ( p,18|) , ( |3.19D and integrating over e we obtain 



Mr 



+ M 



(i) 



.UZ 2 a 3 oo 
9m 2 



L sl 1 + i 



59u L s i 
175u e Li 



3312 / a; L sl 



2401 \u e Lx 



1 

21 



M 4 



here 



+ 9m 2 I 2 , 



l + i 



16cj L s i 384 / L a i 



25tu e Li 245 W e L 



(3.21) 



m 
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A, 2 V V 2 



27rZ 2 a 2 A% a L 1 ' 
Aa, 



1 - f(Za) + 



41 
42 



(3.22) 



Li is defined in ( |3.9|) . The characteristic energy u e encountered in analysis of 
influence of the multiple scattering on the probability of pair photoproduction 
[pm , in gold u; e =10.5 TeV. The amplitudes ( |3.21| ) coincide with the formulas 
( |2.49|) and ([2.571) if we neglect the terms oc u>/u e and (u/u e ) 2 . The terms oc u/u e 
define the real part of the scattering amplitudes while the terms oc {uj/uj e ) 2 are 
the corrections to the imaginary part. 

If the parameter \u s \ > 1 then the value of g c is defined by the equations (see 



±Qq%{qc) 



L, 



L S 3 + - hi 



1 . Aee' 



^s3 — L S 2 + - In L S 2, 
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L s2 = L sl + ~ In — = 2 \n(a s A s ) +1-F 2 (-^ 
Z ujp \ Z 



2f(Za), 



A* = 27rZ 2 a 2 um a 



u Aee' T L 4e£' 
% I L s3 + - In— s- 



(3.23) 



We consider now the region of the strong effect of scattering (\v s \ ^> f ) . We 
restrict to the main terms of the decomposition only, then we can substitute 
the exponential exp(— it) for 1 in integrand in Eqs.( [3.18| ),( |3.19|) . Performing the 
integration we find 



is 2 v s 



FJu H 



is 2 v s ( In 2 — C 



7T 



F 3 {u 8 ) ~ -s 3 u s . (3.24) 



Note, that the next terms of the decomposition can be obtained using the results 
of Appendix A of [pD[ . Substituting these expressions into Eqs. fl3.17|) -( |3.19| ) we 
obtain 



M. 



3am 2 / D L 



s3 



++ 



, 3nZ 2 a*uj 



I f „ 1 7T 

1 2C+ - +i- 

4L s3 V 3 2 



M 4 



2x/2A2 



1/^1 7T 

1 [2C+ - + i- 

4L s3 V 3 2 



(3.25) 



It is seen from these equations that in the case of strong scattering the real and 
imaginary parts of the amplitudes are equal (if we neglect the term oc 1/L s3 in 
M ++ ). Moreover the amplitudes ( |3.25|) don't depend on the electron mass m. In 
place of it we have the value A s . Notice that the asymptotic expansions ( |3.25| ) 
are valid for A C A s (A s > m). In the interval A s 3> A ^> aj 1 we have 



F 2 ~21n(Aa s )-l, F x = 1, L 



s2 



]n^ + l-f(Za) 



(3.26) 



4 Conclusion 



In this paper we considered the influence of a medium on the process of the 
coherent photon scattering illustrated in Fig.l and Fig. 2, where Im M ++ and 
Re M ++ as well as Im M + _ and Re M + _ are given as a function of photon 
energy uj in gold. This influence is due to the multiple scattering of electron and 
positron of the virtual pair on the formation length of the process (see Eqs.( |l.l| ) 
and (|OD) 

U= - (4 1) 

f 2(g s 2 + m 2 + A 2 )' 1 ' ' 

In the region A 2 -C q 2 +m 2 this formation length is independent of A and it value 
is coincide practically with the formation length of pair creation by a photon l c 
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2C]. There is some difference connected with the logarithmic dependence of v. 



2 



value on A 2 : 

,2 



q 2 4nZ 2 a 2 J* dq 2 
Vs\ = —2 = —^ n * l f / -T ( 4 - 2 ) 



J q z 

2 

Q 



where g^ aa , = m 2 +q 2 and = A 2 +a 2 , a is the screening radius of atom ( |1.2|) . 
This defines the weak (logarithmic) dependence of ImM ++ on A in the region 



A < yq 2 + m 2 . This can be seen in Fig.l where the curves 1 and 3 represent 
behavior of Im M ++ for A=0.4435 m and A=0.0387 m respectively. For lower 
value of A the minimal momentum transfer q m in ( |1 ■ 2| ) diminishes thereby the 
interval of contributing the multiple scattering angles increases, so the multiple 
scattering affects the photon scattering amplitude at a lower energy (and smaller 
formation length). Because of this the curve 3 is shifted to the left respect the 
curve 1. Note that the curve 3 (A -1 ~ a s <i) is very similar to the curve 2 in Fig. 2 



of [20] which represents the behavior of the probability of pair photoproduction 
in gold vs photon energy. 

The new property of influence of a medium is the appearance of the real part 
of the coherent photon scattering amplitudes at high energy u. In the region 



to <^ cu e the value of ReM is small accordingly ( |3.21| ). In the asymptotic region 
it) ^> u e we have -Re M = Im M according to ( |3.25|) . This property is seen 
clearly in Figs. 1,2. So the value of -Re M is small at low and very high energies 
of photon. At intermediate energies the value of -Re M have the maximum at 
cu -220 TeV for A =0.4435 m and at uj ~80 TeV for A =0.0387 m. In Fig.2 
the same curves are shown for amplitude M + _. These curves are very similar to 
curves in Fig.l. The curves in both figures are normalized to imaginary part of 
the corresponding amplitude in the absence of the multiple scattering. The ratio 
these imaginary parts one can find from Eqs.( [3.2lD , ( |3.22|) : 

Im Mr _ 1 „ /Aa„\ l\ a, 1 r ' x - ~' 



r = — = —F 1 



Im M ++ I I 



Aa s 



+ 1 



f(Za) + 1 



ln^- - 

A c 2 

(4.3) 

This ratio is r = 0.04435 for A=0.4435 m and r = 0.003018 for A=0.0387 m. 

In numerical calculation of amplitude M ++ we neglect corrections of the order 
1/L S . These corrections are quite small, e.g. for A=0.0387 m they are of the 
order of a few percent. 

We estimate now the integral (over A) cross section of the coherent photon 
scattering at uj u e . From previous section (see Eqs. (|3.25|) - (|3.26|) ) we have for 



module squared of the amplitude at A < A s (see also Introduction) 

\M\ 2 ~ ^ [in + lj , (4.4) 

where value A, is defined in (13.231). It is seen from this formula that in the 



considered region the differential probability of the coherent photon scattering 
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depends weakly (logarithmically) on medium density n a . 



W *( A ) = — v5-"»oA2 = Ai « L ( 4 - 6 ) 



dW ph (A < A t ) = -^\M\ 2 n a dn = i^fdA. (4.5) 

In the interval A ^> A s the influence of the multiple scattering on the coherent 
scattering process is rather weak because the probability of transfer to a medium 
of the momentum > |A| on the formation length If = u/(2A 2 ) appears to be 
small: 

AnZ 2 a 2 to Af 

n = — - 

A 2 a 2A 2 A 4 

In this interval the amplitudes of the coherent photon scattering behave as 1/A 2 
and it doesn't contribute to the integral cross section. So the interval A < A s 
contributes only. Taking into account Eqs. (|4.4j ),( |4.5|) we have for estimate of the 
integral cross section of the coherent photon scattering 

a ph {u » u e ) ~ — '— M A ~ A s 2 ~ -— = 4.7 

lonuj 2 Aj ^27ium, a 
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Figure captions 



Fig.l The amplitude M ++ of the coherent photon scattering in gold under 
influence of the multiple scattering at the different momentum transfer to 
the photon A in terms of the amplitude ImM ++ (|2.57 ) calculated for the 
screened Coulomb potential. 

- Curve 1 is ImM ++ for A =0.4435 m. 

- Curve 2 is ReM ++ for A =0.4435 m. 

- Curve 3 is ImM ++ for A =0.0387 m. 

- Curve 4 is ReM ++ for A =0.0387 m. 

Fig. 2 The same as in Fig.l but for the amplitude M + _ in terms of the 
amplitude ImM + _ ( |2.57| ) calculated for the screened Coulomb potential. 
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Fig. 1 
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Fig. 2 



